Chirality-Dependent Growth Rate of Carbon Nanotubes - A Theoretical Study 



Heiko DumliclB and Stephanie Reich, 

Fachbereich Physik, Freie Universitdt Berlin, 14195 Berlin, Germany 
(Dated: July 20, 2012) 

We consider geometric constraints for the addition of carbon atoms to the rim of a growing nan- 
otube. The growth of a tube proceeds through the conversion of dangling bonds from armchair to 
zigzag and vice versa. We find that the growth rate depends on the rim structure (chirality), the 
energy barriers for dangling bond conversion, and the growth temperature. A calculated chirality 
distribution derived from this minimalistic theory shows surprisingly good agreement with exper- 
iment. Our ideas imply that the chirality distribution of carbon nanotubes can be influenced by 
external parameters. 

PACS numbers: 61.48.De, 61.46.Fg, 81.07.De, 81.10.Aj 



I. INTRODUCTION 

The properties of carbon nanotubes [l| depend 
strongly on their chirality or atomic structure. Most no- 
tably, the metallic and semiconducting character and the 
band gap of a tube change with (n, m) chiral index. One 
of the greatest challenges in nanotube research and ap- 
plication is to control the chirality during the growth. 
This would allow the production of tubes with tailored 
properties without relying on a sorting of bulk samples. 

The growth of a nanotube can conceptually be divided 
into two stages: the nucleation of a cap and the elonga- 
tion of the nucleus into a tube [2-6]. Reich et al. [7] 
showed that the nucleation of the cap fixes the chirality 
of an individual tube as a change in chirality is unlikely 
during the growth phase. Harutyunyan et al. [8] reported 
preferential growth of metallic tubes and claimed the se- 
lection to follow from the shape of the catalytic particles, 
i.e., chirality selection during the nucleation phase. 

The final volume fraction of a given nanotube type does 
not only depend on the nucleation, but also on growth 
speed during elongation. Elongation was mainly studied 
in simplistic models with carbon addition [9, 101. Ding- 
ed al. 11] argued that achiral armchair and zigzag tubes 
grow by introducing kinks when starting a new layer. 
They predicted the armchair kinks to require much less 
energy than zigzag kinks. The growth process, which 
is driven by a monotonous decrease in free energy dur- 
ing elongation, will, therefore, favour armchair tubes. 
Within this line of reasoning chirality selection is inde- 
pendent of external parameters such as catalyst type and 
temperature. 

In this paper, we suggest that the chiral- angle distri- 
bution of carbon nanotubes depends on external param- 
eters. The key is to manipulate the energy difference be- 
tween armchair and zigzag dangling bonds through the 
choice of metal catalyst and growth conditions. We ar- 
rive at this conclusion by looking at the geometry of a 
growing tube, the number and types of places for carbon 
addition. The rim of a nanotube consists of three dif- 
ferent growth sites with a varying energy barrier for the 



addition of carbon atoms. The number and distribution 
of growth sites is a function of chirality. Combining this 
minimalistic geometric approach with calculated energy 
differences for carbon dangling bonds on metals we pre- 
dict a distribution of chiral angles that is in surprisingly 
good agreement with experimental findings. 

This paper is organized as follows. We first show how 
growth proceeds with carbon addition with respect to our 
model, Sec. II A. The essential properties of rims made 
up by hexagons are discussed in Sec. II B. The growth 
factor, which allows us to understand why chiral selectiv- 
ity occurs during the nanotube elongation, is introduced 
in Sec. II C. We then discuss how the chirality distribu- 
tion can be influenced by external parameters in Sec. II 
D. Finally, in Sec. Ill we use our model - derived in Sec. 
II - to obtain two exemplary chirality distributions and 
compare our results to experimentally determined data. 
Section IV summarizes this work. 



II. MODEL 
A. Growth 

Let us first consider schematically the growth of a car- 
bon nanotube. In Fig. [T] we present a 3d wire model 
of a possible growth route of a (5, 5) nanotube. The 
growth proceeds through to addition of C2. The first car- 
bon atom adds endothermically and is followed exother- 
mically by a second carbon atom. The pentagon cre- 
ated in the first step is energetically less favorable than 
hexagons. We, therefore, expect the next carbon atom 
to be added to the pentagon. Also, the creation of more 
and more pentagons would close the tube and terminate 
growth [12|. Alternatively, a carbon dimer is added. 

Going through the series of tubes in Fig. \Tjpi) to d) a 
layer of carbon atoms was grown, which corresponds to 
half a unit cell of the (5, 5) tube. The continuation of the 
process - until growth is terminated - leads to an arm- 
chair carbon nanotube. In the lower panels of Fig-UJa^) to 
d) we introduce a schematic representation of the grow- 
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FIG. 1: Three dimensional (top) and reduced rim (bottom) representations of half layer growth of a (5,5) carbon nanotube. 
The reduced rim representations at the bottom are obtained by unzipping the 3d wire model (top). The "." denotes a growth 
site for C2 addition, the a stands for an armchair and the z for a zigzag dangling bond. The arrow at the left denotes the 
growth direction. We assume root growth, but the picture is turned upside down and a catalyst was omitted for clarity of the 
C2 addition, a) A C2 adds at the rim. The induction of a new layer is accompanied by a barrier [11]. b) 3 C2 additions follow 
without experiencing an energy barrier leading to tube c). The last C2 addition leads to a stable rim or closed layer, d) A half 
layer is grown compared to a). This process continues until the growth is terminated. 



ing rim of a nanotube. The reduced rim representation 
unfolds the rim of the tubes in two dimensions. The "." 
denote growth sites at which the addition of C2 is ener- 
getically favorable because a hexagon is created. The a 
stands for an armchair dangling bond; it consists of one 
of two neighbouring twofold C-C bonded atoms. The z 
stands for zigzag dangling bond; it has two saturated C 
neighbours and is itself twofold C-C bonded. 

The first C2 addition to the rim starts a new layer by 
converting two a into z dangling bonds, which is accom- 
panied by an energy barrier [11[, compare Fig. \flp.) and 
b). The following three C2 additions do not change the 
energy of the rim, since they only move a and z dan- 
gling bonds leading to Fig. \ljp). The final C2 addition 
to Fig. \l]p) yields a finished armchair layer presented in 
d). We argue that the conversion and movement of the 
growth sites and the energetic barriers for the conversion 
determine the chirality-specific growth speed of carbon 
nanotubes. 



B. Rim 

In the following we first discuss the essential properties 
of rims made up by hexagons. All carbon nanotube rims 
consist of armchair a and zigzag- type dangling bonds z. 
The number and nature of the dangling bonds in a rim 
depend mainly on chirality (n,m). The rim that follows 
most closely the chiral circumferential vector of a tube 
has N a = 2m armchair and N z = n — m zigzag dangling 
bonds. For the (5, 5) armchair tube this rim is shown in 
Fig.IUa^. During the growth the total number of dangling 
bonds in the rim remains constant N a + N z = n + m, 
while N a and N z vary. By this condition we include 



TABLE I: Overview of the three growth site types with "A." 
the change in the growth site number, "transition" the change 
of the bond structure and examples for the transitions, 
growth site A. transition example 

aa.aa —1 aa.aa — ► zaaz Fig\l}i)—)>b) 

aa.z/z.aa aa.z — >> zaa./z.aa — »• .aaz ¥\g^))^c) 
z.z +1 z.z — > .aa. Fig\T]p)—>d) 



all reasonable configurations of a growing nanotube and 
exclude obviously unreasonable configurations, e.g., one 
side of the tube being much longer than the other side. 

Following the rims in Figs. \T]p,)-d) from left to right 
leads to the notation of the particular configuration 
of a rim: a) .aa. aa.aa. aa.aa, b) .aa. zaaz. aa.aa, c) 
a.aa.aaz.zaa.a, and d) a. aa.aa. aa.aa. a . A rim thus 
consists of a combination of z.z, aa.z, z.aa and aa.aa 
growth sites. The aa.z and z.aa growth sites are identi- 
cal by symmetry. The aa, za, az, and zz configurations 
do not contain a growth site as C2 addition does not 
add hexagons. Therefore, they do not contribute to the 
tube elongation process. In the starting configuration 
the number of growth sites in a rim of an (n, m) nan- 
otube is N aa . aa = 2m — n, N aa . z = min (m, n — m) with 
"min" the minimum and N ZmZ = 0. Note that for tubes 
with 2m — n < there are only aa.z growth sites and 
zigzag tubes (m = 0) do not contain any growth sites at 
all. The growth of zigzag tubes is supressed in our model 
and needs an intermediate addition of C or C3 to induce 
a growth site for C2 addition, that we do not consider 
here. A C2 addition to the rim will change the type and 
the number of growth sites, see Tab. H 
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C. Growth Factor 

The rim of an (n, m) nanotube with n > m > 
n/2 can be divided into a part with a chiral vector 
(2m — n, 2m — n) that contains aa.aa growth sites and a 
part with a vector (2n — 2m, n — m) that consists exclu- 
sively of aa.z sites. Therefore, all nanotube rims can be 
divided in aa.aa containing rim parts and aa.z contain- 
ing rim parts. During the growth the number of growth 
sites contributed by a rim part containing exclusively 
aa.z sites remains constant. The number of growth sites 
contributed by a rim part containing exclusively aa.aa 
sites, however, changes continously during the growth as 
is best illustrated by the example of an armchair tube. 
Growing a full layer of an armchair rim requires the ad- 
dition of 2 • (2m — n) = n + m carbon dimers. The maxi- 
mum number of aa.aa growth sites 2m — n occurs only at 
half and full armchair layers. The other 2 • (2m — n) — 2 
growth steps have one growth site less[23|. Summing up 
the number of growth sites in each step and dividing by 
the number of C2 additions yields the average growth site 
number 
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Similarily, we find the average growth site number for 
the rim part containing aa.z sites A aa . z = N aa . z . Adding 
the contributions of aa.aa and aa.z rim parts yields the 
average number of growth sites as a function of chiral 
indexes n and m 
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if 2m — n > 
otherwise. 



(2) 



The addition of C2 dimers to the A (n, m) sites will lead 
to a lengthening of the tube with n + m C2 additions 
for a single full layer. If we define the abundance of a 
certain nanotube chirality to depend on the number of 
full carbon layers, we find the growth speed of a tube to 
be proportional to 



T (n, m) 



A (n, m) 
n + m 



(3) 



The growth factor Y (n, m) allows us to understand why 
chiral selectivity occurs during the nanotube elongation 
phase. In the following we will show how we can influ- 
ence the chiral distribution during the elongation of a 
nanotube. 



D. Influence on Chirality Distribution 

The addition of C2 to the different growth sites will 
experience varying energy barriers, as zigzag dangling 



bonds (E z — 2.90 eV) require much more energy than 
armchair dangling bonds (E a — 2.10 eV) in vacuum [l3| . 
The armchair configuration is energetically favorable be- 
cause it consists of two dangling bonds on neighbouring 
C atoms that form a triple bond. To model experimental 
growth conditions we need to consider a metal catalyst 
in most growth scenarios. The energetic difference be- 
tween a and z dangling bonds is reduced by the presence 
of a metal 0,0, El, as carbon- metal bonds are formed. 
However, the difference remains non-zero, as electrons of 
carbon neighbours influence the total bond energy of the 
carbon- metal bonds, rendering a lower in energy than z. 

The energy barrier for the C2 addition to an aa.aa site 
depends on the conversion of aa.aa into zaaz dangling 
bonds (see Tab.QJ. The conversion requires an energy 



E z 



E a 



2E Z - 2E a = 2E a (r - 1) ,(4) 



with E a the energy of an armchair and E z the energy of 
a zigzag dangling bond. With r = E z /E a we denote the 
ratio between the two energies. The total dangling bond 
energies as well as their ratio depend on the catalyst. 
Changing z.z into .aa. we gain A a . In contrast growing 
at an aa.z site will cost no energy; this growth happens 
without an energetic barrier. This energetically different 
behaviour allows to affect the chirality distribution of 
carbon nanotubes through external parameters such as 
the metal catalyst and the growth temperature. 

If the addition to aa.aa sites has a negligible barrier 
(r ~ 1 or A a <C &bT) all growth sites can contribute to 
the growth speed. We combine Eqs. (j3j) and (j2j) to ob- 
tain T. Figure [2]^ shows the growth speed T as area 
size in chiral angle and diameter dependence for diame- 
ters d = 0.675-1.055 nm. The highest V occur for (n,n) 
armchair tubes. A small trend for increasing V exists 
for larger diameter tubes, resulting from the fractional 
term of Eq. (pp), as the comparison of the armchair tubes 
shows. Changing the environment (e.g. another catalyst 
with another r or adjustment of temperature) so that 
A a ^> ksT, the aa.aa growth sites will not contribute 
anymore; Eq. (J3j) yields V = A aCLiZ / (n + m), which leads 
to a different growth speed distribution. The highest T 
now occurs for (n, 77) chiral tubes, see Fig. [2]^. 

For real samples we expect a distribution of growth 
speed T to be between the two limiting cases. The 
thermal energy of nanotube growth is on the order of 
k B T ^ 0.05-0.11 eV 0, E3- A a depends on the cat- 
alyst material, its composition and - less pronounced - 
on the position of the carbon with respect to the metal 
atom. The barriers for metal catalysts are on the order 
of A a & 0-0.12eV for various metals [3, El and thus com- 
parable to the thermal energy. Therefore, the addition 
to the aa.aa site is not suppressed. This agrees with the 
results of Ding et al. [U(, that the barrier for armchair 



kink introduction - which corresponds to C2 addition to 
aa.aa - is negligible. Recently, other materials like Si02 
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FIG. 2: Comparison of T (n, m) for tube diameters d = 
0.675-1.055 nm for a) A a < k B T. b) A a > k B T. The abun- 
dance of metallic/semi- metallic tubes (open circle, red cross) 
decreases compared to semiconducting tubes (full circle, gray 
cross) from a) to b). 



were found to catalyze nanotube growth [16j. Further, 
bimetallic catalysts contain different barriers and may 
be extremely interesting for influencing the chirality dis- 
tribution 



III. DISCUSSION 



d = (0.75 ±0.15) nm and ACCVD samples [18[. The 
agreement between theory and experiment in Fig. [3]^ 
is striking. Our model very well predicts the overall de- 
crease of the number of tubes with increasing chiral angle. 
The strong discrepancies for selected chiralities - e.g. the 
strong luminescence of the (10, 2) tube - is most li kely 
due to a high quantum yield for some nanotubes [22|. 
On the other hand, the nucleation phase might also pre- 
fer certain chiralities [7]. It would be highly desirable 
to establish an unambigious chirality distribution exper- 
imentally to clarify these points. 

Figure EJaj verifies our assumption, that a mixture of 
the T factors derived for A a <C fc#T and A a ^> k B T have 
to be used for real samples as the deviation between the 
theoretical and experimental part of Fig. ^jpt) show. Fig- 
ure ^j)) on the other hand perfectly reproduces the trend 
with only considering the contribution of the T factor for 
A a <C &#T, which is expected to result from the growth 
conditions. We conclude that different growth conditions 
have indeed an influence on the chirality distributions 
which results during the elongation of the nanotubes. 



IV. CONCLUSIONS 

In summary, we suggested how to control the nanotube 
growth and elongation process through the structure of 
the rim. Depending on the tube chirality the rim contains 
three different growth sites. Geometric considerations 
yield the growth factor T, which in turn determines the 
chirality distribution of carbon nanotube samples. We 
showed that chiral selectivity can be obtained through 
a combination of external parameters. Our results will 
be important for the understanding and tailoring of the 
growth process of single-walled carbon nanotubes. 
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Up to now we concentrated on the growth of an ex- 
isting nanotube nucleus. The chirality distribution of 
a sample will also depend on the nucleation phase, i.e., 
whether a particular tube cap is nucleated or not (2o| . 
The diameter of carbon nanotubes is clearly determined 
by the nucleation step [2l[. We now assume a distri- 
bution of chiral indices where (i) the diameter is fixed 
by nucleation and (ii) the chiral angle distribution is 
given by Eq. (|3j) with A a <C k B T. Figure ^fr) com- 
pares the chirality distribution of semiconducting nan- 
otubes with d = (0.93 ± 0.3) nm to the experimental 
distribution in HiPco tubes; Fig. [3]^ is for tubes with 
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